Charge transfer electrostatic model of compositional order in 

perovskite alloys. 



Zhigang Wu and Henry Krakauer 
Department of Physics, College of William and Mary 
Williamsburg, VA 23187 
(February 1, 2008) 



Abstract 



We introduce an electrostatic model including charge transfer, which is shown 
to account for the observed B-site ordering in Pb-based perovskite alloys. The 
model allows charge transfer between A-sites and is a generalization of Bel- 
laiche and Vanderbilt's purely electrostatic model. The large covalency of 
Pb 2+ compared to Ba 2+ is modeled by an environment dependent effective 
A-site charge. Monte Carlo simulations of this model successfully reproduce 
the long range compositional order of both Pb-based and Ba-based complex 
A(BB'b")03 perovskite alloys. The models are also extended to study sys- 
tems with A-site and B-site doping, such as (Na 1 /2La 1 / 2 )(Mg 1 / 3 Nb 2 /3)03, 
(Bai_ 2 La a; )(Mg( 1+a; )/3Nb(2_ a; )/3)03 and (Pbi_ a .La x )(Mg (1+x)/3 Ta ( 2_ x )/3)03. 
General trends are reproduced by purely electrostatic interactions, and charge 
transfer effects indicate that local structural relaxations can tip the balance 
between different B-site orderings in Pb based materials. 

PACS numbers: 64.60.Cn, 81.30.-t, 77.84.-s, 77.84.Dy 
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I. INTRODUCTION 



Complex solid solution perovskite structure based relaxor ferroelectrics such as 
Pb(Zr la; TL E )0 3 (PZT) and (1-ar) Pb(Zn 1/3 Nb 2/ 3)0 3 + x PbTi0 3 (PZN-PT) show great 
promise as piezoelectric transducers. |l] In their simplest form, PZN-PT materials are B- 
site alloys of the form A(BB B )0 3 , with three different B-site cations from group II, IV, 
and V, which exhibit compositionally dependent B-site atomic ordering. For example, at 
1640 °C, when the tetravalent composition x is increased in (l-x) Ba(MgNb)0 3 + x BaZr0 3 
(BMN-BZ) alloy, the following sequence of B-site ordering is observed: [lll]i : 2 order for 
x < 5%; then [111] 1:1 order for 5% < x < 25%; and finally disorder for larger x. ]2[ 
Other Ba-based perovskites, e.g., (l-x) Ba(Mgi /3 Ta 2 / 3 )0 3 + x BaZr0 3 (BMT-BZ), § (1- 
x) Ba(Mgi/ 3 Nb2/3)0 3 + x BaZr0 3 (BMN-BZ), B display a similar sequence of B-site order. 
On the other hand, for Pb-based systems, e.g., (l-x) Pb(Mg!/ 3 Ta 2 / 3 )0 3 + x PbZr0 3 (PMT- 
PZ), [111] 1:2 order is not observed at x = 0; instead, annealing between 1325°C and 1350°C 
results in [111] 1:1 order all the way down to x — 0. |4]|J Other Pb-based perovskites , e.g., 
Pb(Mg 1/3 Nb 2/3 )0 3 (PMN), @-|| display similar B-site ordering. 

Chemical substitutions on the A-site can also strongly afftect B-site ordering. For exam- 
ple, (Nai/ 2 Lai/ 2 )(Mg!/ 3 Nb 2 / 3 )0 3 shows [1 1 l]i : i B-site ordering at all measured temperatures; 
and [001]i : i A-site ordering occurs below ~ 950°C. [Q] In (Bai_ a; La a ;)(Mg( 1+2 .)/ 3 Nb( 2 _ :r )/ 3 )0 3 , 
[lll]i:i B-site ordering is induced and coexists with [111] 1:2 order as low as 5% La doping 
on A-sites. For x > 0.1, [lll]i : i B-site ordering is stablized. Similar B-site ordering se- 



quences are observed for (Bai_ x La x )(Zn(i +x )/ 3 Nb( 2 _ x )/ 3 )0 3 . |T]J On the other hand, anneal- 
ing (Pbi_ x La x )(Mg( 1+a .)/ 3 Ta( 2 _ :r .)/ 3 )0 3 at 1300 °C yields [1 1 l]i : i B-site ordering for all x, and 
the strongest 1:1 reflections are measured at x — 0.5. ]12[ (Pb 1 _ a .La x )(Mg( 1+a ,)/ 3 Nb( 2 _ :r .)/ 3 )0 3 



displays similar B-site ordering sequence. [13 



Bellaiche and Vanderbilt (BV) introduced an electrostatic model that was remarkably 
successful in explaining the B-site ordering of A +2 (BB )0 3 perovskite alloys, using only 
electrostatic interactions between B-site. [|14]] However, the model failed to correctly describe 
the small x behavior of PMT-PZ described above. The BV model has also not been extended 
to study the effects of A-site chemical substitution. 

The different behavior of Ba- and Pb-based alloys may be presumed to arise from the 
greater covalency of Pb compared to Ba. First-principles calculations for Ba based alloys 
show that typical B-site stacking structural energy differences can be as large as 60 kJ/(mol- 
AB0 3 ). ]15| This is an order of magnitude larger than in similar Pb based alloys, [p^-|l7[| and 



reflects the greater covalency of Pb compared to Ba. In this paper, we investigate modeling 
Pb covalency through the inclusion of a local- configuration dependent A-site charge, i.e. 
charge transfer between A-sites. The resulting electrostatic models with and without charge 
transfer are used to study B-site atomic ordering in undoped and La-doped Ba- and Pb-based 
perovskite alloys. 



II. ELECTROSTATIC MODELS WITHOUT CHARGE TRANSFER 

We will be considering supercells of the AB0 3 perovskite structure with possible mixtures 
of different atomic species on the A and B sites. Strain and relaxation of atomic positions 
are neglected in all that follows. The total electrostatic energy is then given by 
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where / indexes the primitive cubic cell within the supercell, r indicates the sublattice 
position in the I th primitive cell (r = {A,B,Oi,02,03}), and e is the electronic dielectric 
constant. The formal oxygen ionic charge Qz,o — — 2 will be used in the following, and the 
A- and B-site charges will be referenced to PbTiCV 

Q lA = 2 + Aq l>A , (2) 

and 

Q ltB = 4 + Aq l:B , (3) 

e.g. Aq t)A = 0, -1, +1 for Ba, Na, and La respectively, and Aq^ B = 0, -2, +5 for Ti, Mg, 
and Nb respectively, give the formal ionic charges for these atoms. BV showed that up to 
a constant, the configurationally averaged total energy depends only on terms quadratic in 
the Aq^ T . The configurational electrostatic energy can then be written as: 
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(4) 



BV introduced this model for A +2 (BB )03 perovskite alloys. Since they only considered 
alloys for which A-sites are occupied by Ba or Pb, where Ag^A = 0, the configurational 
electrostatic energy then depends only on the B-sublattice: 

where a is the cubic lattice constant, and = la. 

This model was remarkably successful in explaining the B-site ordering of A +2 (BB )C>3 
perovskite alloys using Eq. (H). |TJ| Using a B-site-only notation to describe the differ- 



ent classes of alloys simplifies the discussion. Thus, IV X IV 1 _ X denotes a homovalent B-site 
binary alloy having tetravalent B-atoms, e.g., Pb(ZrTi)03; and II( 1 _ a ,)/ 3 IV :r .V2(i_ a ;)/3 indi- 
cates a heterovalent B-site ternary, such as (1-x) Ba(Mg!/ 3 Nb2/3)C>3 + x BaZrCV The 
BV model predicted the following B-site orderings |I4|| : (1) For III( 1 _ a ,)/ 2 IV :r V(i_ 2: )/ 2 and 
II(i-a;)/2lV x VI(i_ a ;)/2 heterovalent ternaries, a rocksalt-type [lll]i : i ordering becomes disor- 
dered with progressively increasing x, consistent with experimental observations. (2) For 
II(i_x)/3lV a; V2(i_x)/3 heterovalent ternaries, the model predicts a richer phase diagram, with 
[lll]i : 2 order for x less than about 5%, [111] 1:1 order for x greater than 5%, followed by 
disorder for large x. This is in good agreement with observations for Ba based systems but 
disagrees for Pb based systems, where [lll]i : i order persists down to x — 0. They proposed a 
possible explanation, []TJJ postulating a small amount of Pb 4+ on the B-sublattice to explain 



the stabilization of the 1:1 phase. As far as we know, there is no evidence of the presence of 
Pb 4+ on the B-sublattice. 
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III. CHARGE TRANSFER ELECTROSTATIC MODEL 



As mentioned, first-principles calculations for Ba based alloys show that typical B-site 
stacking structural energy differences can be as large as 60 kJ/(mol-AB03), [15| which is 



an order of magnitude larger than in the Pb based alloys. ||15|-[T7|1 The energy differences in 
the Ba based compounds can be approximately reproduced in the BV model by choosing 
the value of its one free parameter e = 10 at the experimental volume. [14,H| However, 



the BV model is then unable to account for the much smaller energy differences in the 
Pb based compounds. The greater covalency of Pb compared to Ba permits larger local 
relaxations, reducing the energy differences between different Pb based B-site stackings. 
Covalent interactions are of course absent in a purely electrostatic model. 

We attempt to model the greater Pb covalency by allowing charge transfer between A- 
sites, depending on the local B-site configuration. Our charge transfer electrostatic model 
is closely patterned after similar models that were introduced to model the configuration 
dependence of Coulomb energies for pseudobinary alloys, fi~8^p0| such as (Alo.sGao.5)As. 
These calculations found that Coulomb contributions of calculated mixing enthalpies were 
commensurate with that of lattice mismatch energies. [T^] In addition, the trend of stability 
of (AC)„(BC)„ superlattices predicted by these models was consistent with first-principles 
total-energy calculations for lattice-matched semiconductor superlattices. J2U] 

To avoid confusion in the following, we use an overhead "~" to denote charges that 
depend on the local configuration of a site, while charges without the overhead ~ denote 
ionic charges that depend only on the atom occupying the site. The simplest way to include 
charge transfer in the model is to define a nearest-neighbor configuration-dependent effective 
A-site charge Qi,a'- 

Qi,a = Qa + ^Qi,a = 2 + Aq l:A ■ (6) 

where 

A5U = -A(i £ Ag i>B ). (7) 

i=n. n.B— sites 

The parameter A controls the charge transfer between the A-sites, i. e. the degree of covalency 
on the A-site. Overall charge neutrality is automatically maintained by this expression. With 
this definition the configurationally dependent electrostatic energy now depends on both the 
A- and B-sublattices: 
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This expression is formally similar to Eq. (|]), except that here the charge on an A-site 
depends not only on the atomic species on that site but on the nearest neighbor B-site 
configuration. 

The value of the additional parameter A can be chosen depending on the species occu- 
pying the A-site. Thus for Ba-based alloys A = 0, no charge transfer occurs, and the model 
reduces to Bellaiche and Vanderbilt's model. For Pb-based alloys, the value of A can be fit 
to first-principles calculations of representative supercells of Pb(BB')03. 
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The parameters in our model are chosen as follows, e is chosen to reproduce the energy 
differences of selected Ba-based alloys (we used e = 10, as in Ref. ]T5|). The value of A 
is then fit to selected Pb-based alloy energy differences. Fig. 1 shows structural energy 
differences (x = 0) as a function of A. To obtain B-site stacking structural energy differences 
in the range of 10 kJ/(mol-ABOs) for Pb based alloys energies, as found in first-principles 
calculations, Fig. 1 indicates that Apb should be chosen in the range 0.7 < Apt,. 



IV. MONTE CARLO SIMULATIONS 



Metropolis Monte Carlo simulations |21[ were carried out with periodic boundary con- 
ditions using 6x6x6 supercells, similar to those in BV. To obtain good statistics, typically 
~ 10 7 moves were made at each temperature. Calculations were started at at high tem- 
perature (e.g., 4000 K) with a random configuration, and the temperature was then slowly 
decreased until the acceptance rate becomes very low. To determine the B-site ordering 
the Fourier transform of the charge-charge correlation function r](k) was computed from 
converged runs at a given temperature: WA 



tj(k) = a £ Afl Aq l+l , exp(-ik • l') , (9) 

u' 

where a is a normalization factor, the sum runs over the B-sublattice, and k is the wavevector 
in the Brillouin zone of the unit cubic cell. The charge-charge correlation function r/(k) is 
directly proportional to the ensemble average of the square of the Fourier transform of the 
charge distribution, e.g., a large value of r\ at k = 27r(|, |, |) or k = 27r(|, |, |) corresponds 
to a strong [lll]i : i or [111] 1:2 order respectively. 



V. RESULTS 



MC calculations using electrostatic models with and without charge transfer are reported 
in this section as a function of temperature and concentration x of A- and B-site doping. 
Results for B-site only alloys are presented first. A- and B-site alloy calculations are then 
considered, first using a purely electrostatic models and then including charge transfer. 



A. B-site Alloys With Charge Transfer 



Simulations were performed using Eqs. (|7j,D for A(BB B )0 3 alloys with A = 0.3 and 0.7, 
and setting a = 7.7 a.u. and e = 10. Results for II( 1 _ 2 .)/ 3 IV :c V 2 (i_ x )/3 heterovalent ternaries 
are given in Fig. 2, which shows the calculated charge-charge correlation function r?(k) vs. 
tetravalent composition x at temperatures of T = 1000 K and 1500K and for A = 0.3 and 
0.7 respectively. For A = 0.3, which permits only a small charge transfer between A-sites, 
Bellaiche and Vanderbilt's prediction is reproduced (also for A = 0.0, not shown), i.e. we 
obtain the B-site ordering sequence of Ba-based perovskite alloys described in Section II. 
For A = 0.7, Fig. 2 shows there is no [111] 1:2 ordered phase. Instead, there is a continuous 
transformation from [111] 1:1 order at x = to a disordered state at large x. Increasing the 
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final temperature to T=1500K does not change this ordering sequence, but only decreases 
the value of x at which disorder appears, as expected. 

Moreover, the [lll]i : i structure is similar to the "random-site" model, |2| in that it has 
alternating planes of mixed (Il-V)-planes and pure (V)-planes, as was also found by BV. 
These results are in good agreement with experimental observations of the ordering sequence 
of Pb(Mg 1/3 Nb 2/ 3)0 3 (PMN), p-g (1-x) Pb(Mg 1/3 Ta 2/3 )03 + x PbZr0 3 (PMT-PZ) | and 
Pb(Mg 1/3 Ta 2/3 )0 3 (PMT). § 

The charge transfer model also correctly predicts B-site ordering in the 
ni(i_ 2 )/ 2 IV a; V(i_ a; )/2 and II( 1 _ a .)/ 2 IV x VI(i_ :r )/ 2 alloys. Monte Carlo simulations using Eqs. 

with A = and A = 0.7 show that both Ba-based and Pb-based perovskite alloys 
have the same B-site ordering for III(i_, c )/ 2 rV r V(i_ :r )/ 2 ternaries and 11(i- x )/2NxVIq_- x )/2 
ternaries. Our results reproduce those of BV in these cases. 



B. A- and B-site Alloys Without Charge Transfer 

To study the effects of A-site chemical substitution, we first used a purely electrostatic 
model with formal ionic charges. Calculations were performed for (A^ 2 A 1// 3 2 + )(B^ 3 B 2/ , 5 3 + )0 3 

and (Af1i x A' x 3+ )(B^ +x y 3 B ( ^ x y 3 )03. Since Qz,a is now not a constant, we express the A- 
and B-sublattice charges using Eq. ([|). 

Fig. 3 shows both A- and B-site charge-charge correlation function r/(k) vs. temperature 
for (A}y 2 A 1/ 3 2 + )(B^ 3 B 2/ , 5 3 + )0 3 . The purely electrostatic model, Eq. (|J), can explain the 
observed [111]1:1 B-site ordering of (Na 1/2 La 1/2 )(Mg 1/3 Nb 2/3 )0 3 . g Moreover, the 1:1 B- 
site ordering found in the calculations corresponds to the random site model. As seen in 
Fig. 3, weak [001]i : i A-site ordering is also obtained in the simulations, in agreement with 
measurements on samples annealed below ~ 950°C 0. Perfect [001] m A-site ordering would 
correspond to alternating [001] A-planes of Na and La, but the A-sites in the simulations 
are more disordered (see section VI). In simulations using a 6x6x6 supercell, only B-site 
order was observed. The weak A-site ordering emerged only after increasing the simulation 
region to a larger 12x12x12 supercell. 

Fig. 4 shows 77 (k) vs. tetravalent composition x for (A^ x A :r 3+ )(B^ a ,^ 3 B ( . 2 5 ^,^ 3 )0 3 . The 
symmetry about x = 0.5 is a consequence of our model. Exchanging x with 1 — x in Eq. 
(|), every Aqi T changes sign, and the configurational energy is unchanged. The strongest 
[111]1:1 B-site order occurs at x = 0.5, in agreement with experiment. Experimentally, 
(Bai_ :r La x )(Zn( 1+x )/ 3 Nb( 2 _ a; )/ 3 )0 3 (BLZN) exhibits 1:2 ordering coexisting with 1:1 B-site 



order for x less that ~ 0.05, then switching to 1:1 B-site order for 0.05 < x < 0.6. fil 
(Ba i _ :r La 2: )(Mg( 1+:r )/ 3 Nb( 2 _ a ;)/ 3 )0 3 (BLMN) behaves similarly except that 1:1 B-site order is 
observed over the whole range 0.05 < x < 1. |1(J Fig. 4 shows that the purely electrostatic 
model, Eq. (f|), correctly reproduces the observed B-site ordering, except for the extreme 
La-rich x > 0.95 region, where Eq. (|j) incorrectly predicts 1:2 ordering. However, including 
charge transfer (see below) improves the agreement with experiment in the extreme La-rich 
region (section C below). 

The character of the 1:1 B-site order in Fig. 4 changes at x = 1/2. A detailed 
examination of the equilibrated Monte Carlo structures shows that the 1:1 order re- 
sults from a [111] rocksalt-like alternation of (3' and (3" layers. For x < 0.5 we find 
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[/3 / ]i/ 2 [/3"]i/2 = [(B 2 2 + 2x)/3 )(Bf+ 2x)/3)]i/2[B 5+ ]i/ 2 , i.e. the site is a mixture of B+ 2 and 
B +5 atoms, and the f3" site is solely occupied by B +5 atoms. This is consistent with the 
"random-site" model. The peak in the 1:1 order parameter at x — 0.5 in Fig. 4 corre- 
sponds to [P']x/2 = [B 2+ ]i/2, so the P' and f3" are each occupied by a single cation species, 
and the 1:1 reflections (JToJ exhibit their maximum intensity. However, for 0.5 < x we find 
[/3']i/2[/3"]i/2 = [B 2+ ]i/ 2 [(B 2 2 l x _i)/3)(Bj4 l l2x)/3)]i/2> where now the j3' site is solely occupied by 
B +2 atoms and the j3" site is now a mixture of B +2 and B +5 atoms. Akbas it et al. fI^ , |lT)|Jl"2|| 
proposed this structure based on their observations. 



C. A- and B-site Alloys With Charge Transfer 

We next consider whether the inclusion of charge transfer may improve the agreement 
with experiment for the A- and B-site alloys considered in the previous section. The elec- 
tronegativity and electron affinity of an atom is an indicator of the degree of covalency. By 
this measure, La is more covalent than Ba or Na, and Nb has higher covalency than Ta. 

There is a difficulty in generalizing our charge transfer model to systems with more 
than one atomic species on the A- and B-sites, such as (Bai_ x La a; )(Mg(i +2: )/3Nb( 2 _ :r )/3)03. 
Considering charge transfer only onto La ions and using Eq. (|7|), charge neutrality may be 
violated except for x — 1. For BLMN, the simplest way to restore charge neutrality is to 
permit charge transfer onto both La and Ba ions, using an average charge trasfer parameter 
Aa(x) proportional to the La composition x, i.e., Xa(x) = xAa(1), so that there is no charge 
transfer at x = (BMN) and the largest charge transfer occurs for at x — 1 (LMN). 

In BLMN alloys, charge transfer onto B-sites may also be considered. Mg is pure ionic, 
but Nb and Ta are more polarizable. To ensure charge neutrality we may define Ab similar 
to Aa: 

Ab(x) = ^ T ^A A (1) , (10) 

which varies linearly with B +5 concentration. The configurationally dependent charges on 
the A-sites and B-sites are then defined as: 

Aq lA = -X A (x)(^ Yl A ^b), (11) 

i=n.n.B— sites 



Agi iB = -A B (x)(- Yl A ^a), (12) 

i=n. n. A— sites 

where 

Aqi,r = Qi,r - q T , (13) 

where, r = {A,B}, Q l>T is the formal charge, and <2a and ge are the mean charge of A-sites 
and B-sites respectively. 

We first consider BLMN. To test this model, we set Aa(1) = 0.7, since LMN is observed 
to have [111] i : i B-site ordering. Fig. 5 shows 7/(k) as a function of La composition x for 
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BLMN. The main effect of including charge transfer is that 1:1 order now appears in the 
La-rich region. There is essentially no difference in the B-site ordering with and without 
charge transfer onto the B-sites. Both display [lll]i : 2 order at small x (x < 3%), then 
[lll]i;2 and [1 1 1] i:i order coexist at 3% < x < 5%, and finally [lll]i : i order for all x > 5%. 
However, the model introduces features not seen in the measurements of BLMN. [|1(]] The 
peak of [lll]i:i order shifts from x = 0.5 to x ~ 0.35, and there is a dip in f](k) for [111] i : i 
order at x ~ 0.8. 

The replacement of Ba by Pb in the La-doped perovskites (Pbi- X L& X ) (B/^wgB^J^ , 3 ) 
O3 is now considered. For PLMN, 1:1 B-site ordering is observed for La concentrations, 
increasing continuously up to x = 0.5 and the decreasing for x > 0.5. Again, this trend 
is consistent with Fig. 4, except that the purely electrostatic model results in 1:2 order in 
the extreme Pb and La rich regions. |12[ Several choices for the charge transfer parameter 
were considered. Assuming first that La is much less covalent than Pb, the charge transfer 
constant is chosen to be proportional to the Pb concentration, Aa(x) = (1 — x)\a{1), so that 
there is no charge transfer at x = 1 (LMN) and the largest charge transfer occurs for at x = 
(PMN). The resulting B-site ordering would be similar to Fig. 5, but with x mapped to 1 — x. 
Compared to Fig. 4, this improves the agreement with experiment for the Pb-rich alloys, 
yielding 1:1 B-site order. However, for La concentrations greater than about 95% the model 
does not yield the observed 1:1 ordering. [12] On the other hand assuming that La and Pb 
are similarly covalent, Aa is set to a constant, Aa = 0.7 as in PMN. The results are shown in 
Fig. 6. This model agrees with experiment in the Pb-rich and La-rich regions, reproducing 
the observed 1:1 B-site ordering, but is disordered for intermediate concentrations. Including 
charge transfer onto the B-sites with Xb(x) = 0.7 (not shown) yields similar results to 
those with charge transfer on A-sites only. 



VI. DISCUSSION 

First-principles calculations show that the energy differences between different B-site 
stackings of Pb-based |T5HT7| alloys are greatly reduced compared to those of Ba-based 
fl5| alloys. For example, the energy differences between [100] 1:2, [110] 1:2, [111] 1:2 B-site 
stackings are reduced from ~ 60 kJ/(mol-AB03) in Ba-based alloys to ~ 5 kJ/(mol-AB03) 
in Pb-based alloys. [|T5| The energy reduction can be attributed to the greater covalency of 



Pb compared to Ba, which results in larger local relaxation energies in Pb-based compared 
to Ba-based alloys. As shown in Fig. 1, there is a dramatic decrease in structural energy 
differences when charge transfer between A-sites is included. 

The magnitude of local relaxation energies due to Pb's greater covalency but not includ- 
ing charge transfer between A-sites can be estimated as follows. The ferroelectric PbTiOs 
double-well has a depth of ~ 10 kJ/(mol-ABOs) not including strain, which is much larger 
than the well depth in BaTiOa, consistent with the more covalent nature of Pb. p2[ However 
10 kJ/(mol-AB03) is still not enough to account for the reduction from the 60 kJ/(mol- 
ABO3) B-site stacking energy differences in Ba-based alloys. In heterovalent B-site alloys, 
however, charge transfer between A-sites is another possible local relaxation mechanism. 

This is especially true for small x in B-site ternaries like II(i- a; )/3iV a: V2(i_. c )/3. Since the 
A-site has a nearest neighbor shell of eight B-sites, each of which is occupied by a Ag i B = — 2 
or Aq^B = +1 io n (for x — 0), the sum of these eight Aq^b will necessarily be non-zero in 
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the x = limit. By contrast, in III(i- a! )/2rV'xV(i_ a .)/2 and ll(i- x )/2^x^-(i-x)/2 alloys this 
"non-neutrality" condition is not forced. Our model simulates the effective charge transfer 
between A-sites that results from a combination of oxygen and Pb atomic relaxations and 
direct electronic charge transfer. The purely ionic character of Ba inhibits this, but the 
greater covalency of Pb makes this an effective mechanism to increase the local relaxation 
energy. 

Local relaxation in Pb based alloys plays an important role in compensating the purely 
electrostatic energy differences arising from different B-site stackings. The electrostatic 
charge transfer model in Section III correctly describes the observed B-site ordering in PMN 
type B-site ll(i- x )/3^Vx^2{i-x)/3 heterovalent ternaries. As shown by Fig. 2, introducing 
charge transfer changes the small x ordering from [111] 1:2 to the [111] 1:1 ordering seen 
experimentally. 

While improving agreement with experiment for II-IV-V B-site alloys, the inclusion of 
charge transfer does not alter the correct B-site ordering predictions of the BV model for 
the III-IV-V and II-IV-VI alloys. For these alloys, the simulations predict a rocksalt-like 
structure where the A-site has on average a "neutral" nearest neighbor shell of B-sites. In 
this case, the average Aqa will be small for any value of A. Direct Monte Carlo calculations 
using our charge transfer model predict the same B-site ordering as the purely electrostatic 
model of Bellaiche and Vanderbilt. 

A detailed analysis of our results supports the random-site model of the [111] 1:1 ordered 
regions. We do not find any evidence for the space-charge model. In the space-charge model, 
as discussed in Ref. [|J, the 1:1 ordered regions correspond to a rock-salt like [111] stacking 
of alternating II:V planes, which is locally charge-imbalanced. In the space-charge scenario, 
Coulomb repulsion would limit the size of these charge-imbalanced regions to nanoscale size. 
Overall charge neutrality is restored by embedding the nanoscale 1 : 1 regions in a disordered 
(+5 ion)-rich matrix. The space-charge model was long favored, because it could explain the 
failure of high temperature annealing to increase the size of these 1:1 regions. Instead, the 
recent successful coarsening of these 1:1 regions 0] would seem to rule out charge-imbalanced 
1:1 regions. The random-site model is also supported by recent direct Z-contrast imaging in 
PMN. @ 

Indeed, our calculations show that the [111] 1:1 ordering corresponds to alternating (II- 
V):V planes, i.e. planes with a random mixture of (II- V) cations alternating with pure V 
planes. Bellaiche and Vanderbilt also observed this structure in their 1:1 ordered regions. 
This stacking is consistent with the random-site model. In materials that are also A-site 
alloys, our calculations reveal interesting variants of the random-site model, as discussed 
below. 

The results of A- and B-site alloy B-site ordering presented in the previous section show 
that the purely electrostatic model, Eq. (£|), can account for many of the observations. 
For (Nai/ 2 La 1 /2)(Mg 1 / 3 Nb 2 /3)03 Dupont et al. || argued that [111] 1:2 B-site order in is 
incompatible with the different valences of the A-site cations. In the [111] 1:2 structure there 
are two types of environments for the A-site. In the first, two A-sites are located between 
a B 5+ and a B 2+ layer and have 5 B 5+ and 3 B 2+ nearest neighbors. The formal charge on 
these subcells is [(5(+5) + 3(+2))/8 — 6] = —2.125. The second A-site is located between 
two B 5+ layers and has 6 B 5+ and 2 B 2+ nearest neighbors, with a subcell formal charge of 
-1.75. With a 1:1 mixture of Na +1 and La +3 on the A-sites, each cation must occupy both 
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types of A-site, leading to local charge imbalances of -1.125 and -0.75 for Na and +0.875 
and +1.25 for La on the two respective type of A-sites. These are much larger than the 
charge imbalances when all A-site are occupied by divalent cations, namely -0.125 and +0.25, 
respectively. By contrast, in the random-site 1:1 B-site order structure, all the A-sites are 
equivalent and have 4 B 5+ and 4 (B 2 1g Bff 3 ) neighbors, with a formal subcell charge of -2. 
This leads to local charge imbalances of -1 for Na and +1 for La. Dupont et al. noted that 
this reduction could be a factor in promoting the [111] 1: 1 B-site ordered structure. The MC 
results in Fig. 3, using the electrostatic-only model, confirm that 1:1 order can arise from 
purely electrostatic interactions. 

As mentioned, the weak A-site ordering was not seen in simulations using a 6x6x6 
supercell, but only in the larger 12x12x12 supercell shown. The large error bars in Fig. 3 
reflect the difficulty in equilibrating this system at the lower temperatures in the Monte Carlo 
runs. This is probably due to the presence of many structural states with nearly degenerate 
energies. Perfect [001] i : i A-site ordering would correspond to alternating [001] A-planes of 
Na and La, but this is not seen in the simulations and is reflected in the small value of the 
order parameter. To test if perfect [001] \ : \ A-site ordering is a global energetic minimum, 
we constrained the A-sites to have this ordering and allowed only the B-site configuration to 
change. This resulted in much weaker [111] 1:1 B-site order than seen in Fig. 3 and indicates 
the complex nature of the A-site order. 

In BLMN, Fig. 4 shows that purely electrostatic interactions alone can account for 
stabilization of 1:2 B-site ordering for small La concentrations, switching to 1:1 ordering 



for larger concentrations. The experimentally observed 1:1 reflections |T0[ also have a peak 
at x = 0.5, which is consistent with the random-site model. For x < 0.5 we find the j3' 
site is occupied by a random mixture of B +2 and B +5 atoms, while for x > 0.5 the j3" 
site is occupied by a random mixture. At x = 0.5, the (3' and (3" are each occupied by a 
single cation species, yielding a peak in the 1:1 reflections. Experimentally, the 1:1 ordering 
persists up to the extreme La rich region |Tt]]. However, the intensity of the 1:1 ordering 
reflections decreases between 0.5 < x < 1.0, as reproduced by Fig. 4. Akbas and Davies 
have argued that the 1:2 structure structure is stabilized by anion displacements. In BMN, 
1:2 MNN ordering into separate B +2 and B +5 layers permits adjacent anion layers to undergo 
a concerted displacement towards the smaller B +5 ion. This requires Nb cations to have an 
asymmetric anion coordination. For x = 1 (LMN), 1:2 MMN ordering would lead to an 



asymmetric environment for Mg, which is rarely observed. |10[ Akbas and Davies argued 
that this stabilizes the 1:1 ordered phase for x = 1. The present calculations show that 
purely electrostatic interactions alone stabilize the 1:2 ordering near x = and x = 1, since 
atomic displacements are not included. The results in Fig. 5 and 6 show that 1:1 order at 
x = 1 can be stabilized by the inclusion of charge transfer, which serves to reduce the energy 
differences between different B-site orderings, as shown in Fig. 1. For BLMN and PLMN, 
we interpret this charge transfer as modeling larger local structural relaxations for La than 
for Ba, tipping the balance slightly in favor of 1:1 ordering. 



VII. CONCLUSIONS 

We have used electrostatic models with and without charge transfer to study B-site 
ordering in materials with alloying on both the A- and B-sites. For Pb 2+ and possibly other 
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covalent ions, the effective charge on the A-site can depend on the local B-site configuration. 
The inclusion of charge transfer into the electrostatic model of Bellaiche and Vanderbilt is 
shown to reproduce the observed B-site ordering in PMN type ternary perovskite alloys. For 
A- and B-site alloys, general trends are reproduced by purely electrostatic interactions. In 
(Nai/ 2 Lai/ 2 )(Mgi/ 3 Nb2/3)03, electrostatic effects alone produce B-site and A-site ordering 
observed in experiments. Charge transfer effects in simulations for Pb-based A- and B-site 
alloys indicate that local structural relaxations can tip the balance between different B-site 
orderings in some cases. With the inclusion of charge transfer in these cases, we conclude 
that the electrostatic interaction is the leading mechanism responsible for the B-site ordering 
in perovskite alloys. 
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FIGURES 




Fig. 2 Wu and Krakauer 
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FIG. 2. Monte Carlo simulations of the long-range order parameter ??(k) vs. tetravalent atomic 
composition x using the charge transfer model for II( 1 _ a .)/ 3 IV x V2(i_ a; )/3 ternaries, with 6x6x6 
supercells, a = 7.7 a.u., and e = 10. 
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Fig. 3 Wu and Krakauer 
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FIG. 3. Monte Carlo simulations of the long-range order parameters r/(k) vs. temperature T 
using BV model for (A}/ 2 A 1/ , 3 2 + )(B^ 3 B 2/ , 5 3 + )03 perovskites, with 12x12x12 supercells, a = 7.5 a.u., 
and e = 10. 



Fig. 4 Wu and Krakauer 




FIG. 4. Monte Carlo simulations of the long-range order parameter r/(k) vs. A' x 3+ composi- 
tion x using BV model for (Af^ x A x 3+ )(B^ +x y 3 B^ x y 3 )03 perovskites, with 6x6x6 supercells, 
a = 7.5 a.u., and e = 10, at temperature T = 1000 K. 
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Fig. 5 Wu and Krakauer 




FIG. 5. Monte Carlo simulations of the long-range order parameter 77 (k) vs. La composition x 
using charge transfer model for (Bai-zLa^B?^ n / 3 B , 3 )03 perovskites, with 6x6x6 super- 
cells, a = 7.5 a.u., and e = 10, at temperature T = 1000 K, X\(x) = 0.7x, Xb(x) = 0.7(2 — x)/3. 
The solid symbols denote charge transfer only on A-sites, and the unfilled symbols denote charge 
transfer on both A-sites and B-sites. 



Fig. 6 Wu and Krakauer 



• • [111]1:1 

g — - o [111]1:2 




0.0 0.2 0.4 0.6 0.8 1.0 



X 

FIG. 6. Monte Carlo simulations of the long-range order parameter r?(k) vs. La composition x 
using charge transfer model only on A-sites for (Pbi_ x La x )(B/^ s , 3 B, 2 ^w 3 )03 perovskites, with 
6x6x6 supercells, a = 7.5 a.u., and e = 10, at temperature T = 1000K. Aa = 0.7. 
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